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The recent developments on PYQUEN, HYDJET and HYDJET++ event generators are presented. 
The partonic energy loss model PYQUEN is implemented as the modification of the "standard" 
jet event obtained with the generator of hadron-hadron interactions PYTHIA. HYDJET and HY- 
DJET++ are the Monte-Carlo event generators for the simulation of relativistic heavy ion AA 
collisions considered as a superposition of the soft, hydro-type state and the hard, multi-parton 
fragmentation. HYDJET++ model is the development and continuation of HYDJET generator, 
and it includes more detailed treatment of the "thermal" hadronic state generated on the chemi- 
cal and thermal freeze-out hypersurfaces (represented by the parameterization of relativistic hy- 
drodynamics with preset freeze-out conditions), collective flow effects and decays of hadronic 
resonances. Some applications of above models to novel jet quenching observables are discussed. 
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1. Introduction 

The experimental and phenomenological study of multi-particle production in relativistic heavy 
ion collisions is expected to provide valuable information on the dynamical behaviour of strongly- 
interacting matter in the form of a quark-gluon plasma (QGP), as predicted by lattice Quantum 
Chromodynamics (QCD) calculations. Ongoing and future experimental studies in a wide range 
of heavy ion beam energies require the development of new Monte-Carlo (MC) event generators 
and improvement of existing ones. Especially for experiments at the CERN Large Hadron Col- 
lider (LHC), because of very high parton and hadron multiplicities, one needs fast (but realistic) 
MC tools for heavy ion event simulations Q]. The main advantage of MC technique for 

the simulation of high-multiplicity hadroproduction is that it allows a visual comparison of the- 
ory and data, including if necessary the detailed detector acceptances, resolutions and responses. 
The realistic MC event generator has to include maximum possible number of observed physical 
effects, which are important for the determination of event topologies: from bulk properties of 
soft hadroproduction (domain of low transverse momenta pj < lGeV/c), such as collective flows, 
to hard multi-parton production in hot and dense QCD-matter, which reveals itself in spectra of 
high- pj particles and hadronic jets. Moreover, the role of hard and semi-hard particle produc- 
tion at LHC can be significant even for bulk properties of created matter, and hard probes of QGP 
became clearly observable in various new channels [Q]. In most of the available MC heavy ion 
event generators, the simultaneous treatment of collective flow effects for soft hadroproduction and 
hard multi-parton in-medium production (medium-induced partonic rescattering and energy loss, 
so called "jet quenching") is lacking. Thus, in order to analyze existing data on low and high-pr 
hadron production, to test the sensitivity of physical observables at the upcoming LHC experiments 
(and other future heavy ion facilities) to the QGP formation, and study the experimental capabil- 
ities of constructed detectors, the development of adequate and fast MC models for simultaneous 
collective flow and jet quenching simulations is necessary. HYDJET and HYDJET++ event gener- 
ators includes detailed treatment of hard multi-parton production with taking into account known 
medium effects (such as jet quenching and nuclear shadowing), as well as soft hadroproduction. 



2. PYQUEN event generator 

The event generator for medium-modified nucleon-nucleon collisions PYQUEN ^ is the 
modification of the "standard" jet event obtained with the generator of hadron-hadron interac- 
tions PYTHIA The details of PYQUEN partonic energy loss model can be found in [[7|]. The 
event-by-event simulation procedure in PYQUEN includes the following steps. 

1. Generation of initial partonic state production by calling PYTHIA (parton fragmentation be- 
ing switched off). The spatial vertex of a jet production J(x = rcos y,y = rsin y) is generated 
according to the distribution 

W T t (n)-W i (21) 



dyrdr T^b) 
where b is the impact parameter of a heavy ion collision, Ta and Taa are the standard nuclear 
thickness and nuclear overlap functions, ri^(b,r,Y) are the distances between the nucleus 
centers and the vertex J, and r is the distance from the nuclear collision axis to vertex J. 
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2. Then initial PYTHIA-produced partonic state is modified by medium-induced radiative and 
collisional energy loss. For each hard partem the following rescattering scheme is applied in 
the course of a proper time T. 

• Calculation of the scattering cross section cj(t ; ) = Jdt da/dt and generation of the trans- 
verse momentum transfer t(x{) in the i-th scattering in the high-momentum trasfer limit. 

• Generation of the displacement between the z'-th and (/+ l)-th scatterings, = (t,- + i — T,): 

k 

dP f 

— =X-\z i+1 )exp(- jX-\T i + s)ds), X-\x) = a(x)p(x) , (2.2) 

! o 

and calculation of the corresponding transverse distance, lipj/E. Here A = I /(op) is in- 
medium mean free path, p <=c r 3 is the medium density at the temperature T, o is the integral 
cross section for parton interaction in the medium. 

• Generation of the energy of a radiated in the z'-th scattering gluon, ft), = AE Ta d,i, according 
to BDMS radiation spectrum [M |9|] and its emission angle 0, relative to the parent parton 
determined according to the selected by user parameterization. 

• Calculation of the collisional energy loss in the z'-th scattering 

A£coi,; = ^- , (2.3) 
zrao 

where energy of "thermal" medium parton mo is generated according to the isotropic Boltz- 
mann distribution at the temperature T(Xj). 

• Reducing the parton energy by collisional and radiative loss per each z'-th scattering, 

A£ toM = AE cohi + AE mdJ , (2.4) 

and changing the parton momentum direction by adding the transverse momentum kick due 
to elastic scattering i, 

K = (E--^f-(p---^-±) 2 -ml. (2.5) 
2m 0i p2m 0i 2p 1 

• Going to the next rescattering, or halting the rescattering if one of the following two con- 
ditions is fulfilled: (a) the parton escapes the hot QGP zone, i.e. the temperature in the next 
point r(T i+ i,r !+ i,T], + i) becomes less than T c = 200 MeV; or (b) the parton loses so much 
energy that its transverse momentum pj(Xi+i) drops below the average transverse momen- 
tum of the "thermal" constituents of the medium, 2T(x i+ i). In latter case, such a parton is 
considered to be "thermalized" and its momentum in the rest frame of the quark-gluon fluid 
is generated from the random "thermal" distribution, dN/d 3 p oc exp (-E/T), boosted to the 
center-of-mass of the nucleus-nucleus collision. 

3. At the end of each NN collision, adding new (in-medium emitted) gluons to the PYTHIA par- 
ton list and rearrangement of partons to update string formation with the subroutine PYJOIN 
are performed. An additional gluon is included in the same string as its "parent", and colour 
connections of such gluons are re-ordered by their z-coordinates along the string. Then final 
hadrons are formed by calling standard PYTHIA fragmentation routine PYEXEC. 
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3. HYDJET and HYDJET++ event generators 

HYDJET [@, [KJ and HYDJET++ |JTTJ, ||] are the Monte-Carlo event generators for the simu- 
lation of relativistic heavy ion AA collisions considered as a independent superposition of the soft, 
hydro-type component and the hard, multi-parton fragmentation. The hard parts of HYDJET++ 
and HYDJET are identical, but soft parts are different. The treatment of soft component in HY- 
DJET is rather oversimplified. It considers pions, kaons and protons/neutrons only, the hadron 
spectrum being generated as the convolution of thermal motion and collective flow [|7|]. The soft 
part of HYDJET++ is the "thermal" hadronic state generated on the chemical and thermal freeze- 
out hypersurfaces represented by the parameterization of relativistic hydrodynamics with preset 
freeze-out conditions(the adapted C++ code FAST MC [|13|, 14]). It includes longitudinal, radial 



and elliptic flow effects and decays of hadronic resonances. The HYDJET was written in Fortran. 
The main program of HYDJET++ was written in C++ under the ROOT environment (but there 
is the Fortran-written hard part included in the generator structure as a separate directory). 

Before any event generation, HYDJET/HYDJET++ run starts from PYTHIA initialization 
at given c.m.s. energy per nucleon pair yfs, and then it calculates the total inelastic o™ N {^) 
and hard scattering o^ff L {y/s,p™ m ) NN cross sections with the minimum transverse momentum 
transfer pf m . If the impact parameter b of heavy ion AA collision is not fixed by user, its value b is 
generated in each event between minimum and maximum values in accordance with the differential 
inelastic AA cross section: 

A 2 



a a ™ (b Jj) 
d 2 b { ,V J 



(3.1) 



If the impact parameter b is fixed, then its value in each event is equal to the corresponding input 
parameter. After specification of b for each given event, mean numbers of binary NN sub-collisions 
A^bin and nucleons-participants _/V part are calculated: 

Wn{b,V~s) = TMWofafiVs) , (3.2) 
2n ~ 

iVfa, Vj) =J dx > / J rdr Mn) [1 -exp{o^ N (Vs)TA(r 2 )}] ■ (3.3) 


The next step is the simulation of properly particle production in the event. The soft, hydro-type 
state and hard, multi-parton state are simulated independently. When the generation of soft and hard 
components in each event at given b is completed, the event record (information about coordinates 
and momenta of primordial particles, decay products of unstable particles and stable particles) is 
formed as the junction of these two independent event outputs. 

3.1 Generation of hard multi-parton state 

1 . Calculation of the number of NN sub-collisions Nj^ producing hard parton-parton scatterings 
of selected type with pj > pf m , according to the binomial distribution around its mean 
value (without nuclear shadowing correction yet, 5=1). For this purpose, each from A^m 
sub-collisions is treated by the comparison of random number ^, generated uniformly in the 
interval [0,1] with the probability P^d = ^NN d / a NN produce the hard process. The i-th 
sub-collision is accepted if < Phard, and is rejected in the opposite case. 
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2. Selecting the type of hard NN sub-collision (pp, np or nn) in accordance with the phenomeno- 
logical formula for number of protons Z in the stable nucleus A, Z = A/ (1.98 + 0.015A 2 / 3 ). 
For this purpose, every from Nj et "successful" sub-collisions is treated by the comparison 
of two random numbers and generated uniformly in the interval [0, 1] with the prob- 
ability Z/A. The proton-proton sub-collision is selected if < Z/A, neutron-neutron 
sub-collision — if , > Z/A, and proton-neutron sub-collision — in other cases. 

3. Generation of multi-parton production in jtyj et hard NN sub-collisions by calling PYQUEN 
Afjet times (see the previous section). 

4. If nuclear shadowing is switched on and beam ions are Pb, Au, Pd or Ca, each hard NN 
sub-collision is treated by the comparison of random number 5 generated uniformly in the 
interval [0, 1] with the shadowing factor S taken from the available parameterization Jl6|]. It 
is determined by the type of initially scattered hard partons, momentum fractions taken by 
the partons at the initial hard interaction x\£, the square of transverse momentum transfer 
in the hard scattering Q 2 , and the transverse position of jet production vertex relative to the 
the nucleus centers r^. The given sub-collision is accepted if < S, and is rejected in the 
opposite case. 

5. Formation of hadrons for each "accepted" hard NN sub-collision with PYTHIA (parton frag- 
mentation being switched on), and final junction of jtyj et sub-events to common array using 
standard PYTHIA event output format. 



3.2 Generation of soft "thermal" state 

The details of physics models used to generate thermal hadronic states in HYDJET and HY- 
DJET++ can be found in [J7|] and in [11] respectively. The following MC simulation procedure is 
applied to generate the "soft" thermal component in HYDJET++ (in HYDJET only step (4) for 
pions, kaons and protons/neutrons is used). 



1. Initialization of the chemical freeze-out parameters. It includes the calculation of particle 
number densities. So far, only the stable hadrons and resonances consisting of u, d, s quarks 
are taken into account from the SHARE particle data table [|l7|]. 

2. Initialization of the thermal freeze-out parameters (if T th < T ch ). It includes the calculation 
of chemical potentials and particle number densities. 

3. Calculation of the effective volume of hadron emission region V e ff(b), the fireball transverse 
radius R/(b), the freeze-out proper time X/(b) and emission duration Ax/(b), and the mean 
multiplicity of each particle species. Then the multiplicity is generated around its mean value 
according to the Poisson distribution. 

4. For each hadron the following procedure to generate its four-momentum is applied. 

• Generation of four-coordinates of a hadron in the fireball rest frame 

= {icoshT], rcos0, rsin0, TsinhT]} on each freeze-out hypersurface segment x(r) for 
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the element d 3 o ll u tl = d 3 a^ = n* (r) | 1 — (dx/dr) 2 \ l l 2 x(r)d 2 rdr\, assuming n*Q and X func- 
tions of r (i.e., independent of v\ , <j>). It includes sampling uniformly distributed <p in the 
interval [0, 2n], generating v\ according to the Gaussian distribution exp(— T] 2 /2T]^ ax ) and 
r in the interval [0,Rf(b)]) using a 100% efficient procedure similar to the ROOT routine 
GetRandom(). 

• Calculation of the corresponding collective flow four-velocities. 

• Generation of the three-momentum of a hadron in the fluid element rest frame 

p* {sin d*cos(j)*, sin d*sin(l)*, cos0*} according to the equilibrium distribution function 
f P (p°* ; T, Hi) p* 2 dp*d cos 6* d(j)* by sampling uniformly distributed cos0* in the interval 
[—1,1] and </>* in the interval [0,27i], and generating p* using a 100% efficient procedure 
(similar to ROOT routine GetRandom()). 

• The standard von Neumann rejection/acceptance procedure to account the difference be- 
tween the true probability W* j d 3 ad 3 p*/p°* and the probability n°*jf q (p°* ; T, \ii)d 2 fdr\d 3 p* 
corresponding to the previous simulation steps. For this purpose, the residual weight is cal- 



culated [13]: 



W*d 3 o ( n*p* \ 

W™ = - = x 1 - 4-V . (3.4) 

Then the simulated hadron four-coordinate and four-momentum is treated by the comparison 
of W ( - res with the random number ^ generated uniformly in the interval [0,max(W ( res )]. The 
j-th hadron is accepted if ^ < W, res , and is rejected in the opposite case (then the generation 
of its four-coordinate and four-momentum is tried again). 

• Boost of the hadron four-momentum in the center mass frame of the event using the velocity 
field v(x), that is, 

p ° = Y (p°*+vp *), P = P* + yi} + 7)-\p°* +p°)v . (3.5) 

Note that a high generation speed for this algorithm is achieved because of almost 100% 
generation efficiency due to nearly uniform residual weights W[ ss . 

5. Formation of final hadrons from the two- and three-body decays of resonances with the ran- 
dom choice of decay channels according to the branching ratios taken from SHARE particle 
data files. The two- and three- particle "decayer" programs were developed by N.S. Amelin 
for FAST MC generator [|3[ and have been implemented in HYDJET++. 

4. Physics applications of PYQUEN and HYD JET/HYD JET++ 

It was demonstrated in [^] that PYQUEN and HYDJET models are capable of reproducing the 
main features of jet quenching pattern in Au+Au collisions at RHIC (high-/?r hadron spectra and 
the suppression of azimuthal back-to-back correlations). HYDJET++ describes well even wider 
range of various hadronic observables measured in heavy ion collisions at RHIC for different cen- 
trality sets and kinematic ranges: ratio of hadron yields, pseudorapidity and transverse momentum 
spectra (for both low- and high-p^ domains), elliptic flow coefficients V2(pr), femtoscopic corre- 



lations [11]. Thus the above generators can be applied for the various simulation studies at LHC 



energies. The set of PYQUEN/HYDJET applications to novel jet quenching observables in heavy 



6 



Recent developments on HYDJET/PYQUEN event generators and novel jet quenching observables at LHC 
I.P. Lokhtin 



ion collisions at the LHC was summarized in [Qj. Let us enumerate the main predicted observables 
(all numbers were obtained for central Pb+Pb collisions at y/s = 5.5A TeV and hadronic jets of a 
cone size R = 0.5 and transverse energy Ej > 100 GeV). 

Nuclear modification factor for jets. The jet nuclear modification factor is defined as a ratio of 
jet yields in AA and pp collisions normalized on the number of binary nucleon-nucleon collisions. 
The predicted by PYQUEN/HYDJET jet suppression factor (due to partial gluon bremsstrahlung 
out of jet cone and collisional loss) is about 2 and almost independent on jet energy. It is clear that 
the measured jet nuclear modification factor will be very sensitive to the fraction of partonic energy 
loss carried out of the jet cone. 

Medium-modified jet fragmentation function. The "jet fragmentation function" (JFF), D(z), 
is defined as the probability for a given product of the jet fragmentation to carry a fraction z of 
the jet transverse energy. Significant medium-modified softening of the JFF (by a factor of ~ 4 
and slightly increasing with z) was obtained with PYQUEN/HYDJET simulations. In addition, 
the "anti-correlation" between softening of the JFF and suppression of the absolute jet rates due 
to partonic energy loss out of jet cone (wide-angle gluon bremsstrahlung and collisional loss) is 
predicted. 

Jets induced by heavy quarks. The possibility to observe the medium effects in the channel with 
the production of B-jet tagged by an energetic muon was also analyzed. A significant softening of 
the b-jet fragmentation function (measured with the energetic muon) due to b-quark energy loss in 
QGP is predicted. 

Azimuthal anisotropy of jet production. The azimuthal anisotropy of particle spectrum is char- 
acterized by the second coefficient of the Fourier expansion of particle azimuthal distribution, el- 
liptic flow coefficient, v%. The non-uniform dependence of medium-induced partonic energy loss in 
non-central heavy ion collisions on the partem azimuthal angle cp (with respect to the reaction plane) 
is mapped onto the final jet spectra. The predicted by HYDJET/PYQUEN elliptic flow coefficients 
for jets v 2 increases almost linearly with the growth of b and becomes a maximum at b ~ L.oRa 
(where is the nucleus radius). After that, the coefficients drop rapidly with increasing b. The 
maximum estimated value is v^ 1 ~ 0.04. 

Pj -imbalance in dimuon tagged jet events. An important probe of medium-induced partonic 
energy loss in ultrarelativistic heavy ion collisions is production of a single jet opposite to a gauge 
boson such as y*/Z° decaying into dileptons. The advantage of such processes is that the mean 
initial transverse momentum of the hard jet equal to the mean initial/final transverse momentum of 
boson, and the energy lost by the parton can be estimated from the observed p T -imbalance between 
the leading particle in a jet and the lepton pair. Despite the fact that the initial distribution is smeared 
and asymmetric due to initial-state gluon radiation, hadronization effects, etc., the predicted by 
HYDJET/PYQUEN additional smearing and the displaced mean and maximum values of the p-\- 
imbalance due to partonic energy loss can be significant. The /^-imbalance between the jU + jU~ pair 
and a leading particle in a jet is directly related to the absolute value of partonic energy loss, and 
almost insensitive to the form of the angular spectrum of the emitted gluons and to the experimental 
jet energy resolution. 
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5. Summary 

Among other heavy ion event generators, H YD JET and HYDJET++ are concentrated on the 
detailed simulation of jet quenching effect based on the partonic energy loss model PYQUEN, and 
also reproducing the main features of nuclear collective dynamics by fast (but realistic) way. The 
final hadron state in HYDJET/HYDJET++ represents the superposition of two independent com- 
ponents: hard multi-parton fragmentation and soft hydro-type part. The hard parts of HYDJET and 
HYDJET++ are identical. The soft part of HYDJET++ contains the important additional features 
as compared with HYDJET (detailed treatment of thermal and chemical freeze-out hypersurfaces, 
resonance decays, etc.) HYDJET describes well the main high-pr observables in heavy ion col- 
lisions at RHIC, while HYDJET++ is capable of reproducing also the bulk event properties (such 
as hadron spectra and ratios, radial and elliptic flow, femtoscopic momentum correlations). Both 
generators are promising tools for the various simulation studies in heavy ion collisions at the LHC. 
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